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The strongly correlated spin-electron system on a diamond chain containing localized Ising spins 
on its nodal lattice sites and mobile electrons on its interstitial sites is exactly solved in a magnetic 
field using the transfer-matrix method. We have investigated in detail all available ground states, 
the magnetization processes, the spin-spin correlation functions around an elementary plaquette, 
fermionic quantum concurrence and spin frustration. It is shown that the fermionic entanglement 
between mobile electrons hopping on interstitial sites and the kinetically-induced spin frustration 
are closely related yet independent phenomena. In the ground state, quantum entanglement only 
appears within a frustrated unsaturated paramagnetic phase, while thermal fluctuations can pro¬ 
mote some degree of quantum entanglement above the non-frustrated ground states with saturated 
paramagnetic or classical ferrimagnetic spin arrangements. 
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I. INTRODUCTION 

The entanglement of physical states unveils the pres¬ 
ence of non-local correlations in quantum spin systems 
that have no classical counterpart M. Recently, the 
theoretical and experimental study of quantum entangle¬ 
ment has been put forward by a possibility of informa¬ 
tion processing at the quantum level Q . The generation 
and manipulation of entangled states are fundamental in 
quantum information processes such as quantum compu¬ 
tation [El , teleportation [bl, 0 , cloning of quantum states 
[H and quantum cryptography [9|. Along this direction, 
condensed matter systems play a central role in the study 
of quantum entanglement, because several quantum de¬ 
vices have been proposed on the basis of solid-state sys¬ 
tems. It has been demonstrated that quantum entangle¬ 
ment can affect the low-temperature behavior of macro¬ 
scopic properties, such as the magnetic susceptibility and 
specific heat [iMl. This fact has raised the interest 
in exploring the relation between quantum ent ang lement 
and thermodynamic macroscopic observables |1^ . 

Within the above scenario, quantum spin chains serve 
as ideal model systems in studies of entanglement signa¬ 
tures [iM3. In particular, it has been shown that quan¬ 
tum entanglement exhibits a characteristic scaling behav¬ 
ior in a vicinity of the quantum critical point present in a 
class of one-dimensional quantum spin systems [TsI . This 
feature has opened a new perspective for the description 
of quantum phase transitions |19l-[^ and resulted in the 
introduction of new quantities to measure the degree of 
quantum entanglement such as the quantum concurrence 

The quantum Heisenberg spin chain is the prototype 
model used to search for entanglement signatures in the 
thermodynamic properties of magnetic systems min- 
1^ . In an isotropic Heisenberg system composed of two 
qubits, the quantum concurrence decreases as a function 


of temperature until it completely vanishes at a char¬ 
acteristic threshold temperature, which depends on the 
applied external magnetic field. Further, the degree of 
quantum entanglement decreases monotonically with in¬ 
creasing the magnetic field for all t emp eratures . On 
the other hand, Star ace et. al. showed that the 

anisotropy and external field can be used to produce and 
control in a two-qubit XY Heisenberg chain the degree 
of quantum entanglement at any finite temperature. 

In the last two decades, quantum diamond spin chains 
have been also subject to extensive investigations. The 
spin-1/2 Heisenberg diamond chain |29l - l^ as well as 
its simplified Ising-Heisenberg version j38l-l^ have been 
largely explored in connection with a possible interplay 
between geometric spin frustration and quantum fluc¬ 
tuations. These systems display a rich variety of un- 
usu al p hysical properties such as magnetization plateaus 
[35I . . doubled peaks in the specific heat and suscep¬ 

tibility [H, [sil and quantum entanglement • In addi¬ 
tion, magnetic properties of A3Cu3(P04)4 with A=Ca,Sr 

IS IS, Cu3(Te03)2Br2 lii, K3CU3AIO2(804)4 S, 

Cu 3(C03)^H)2 (azurite) H and Cu3(0H)5(N03) 
(likasite) are satisfactorily captured by various ver¬ 
sions of the spin-1/2 Heisenberg diamond chain. High- 
held magnetization measurements on the natural miner¬ 
als azurite and likasite have for instance conhrmed an 
existence of one-third [46| and two-thirds [13 magneti¬ 
zation plateau besides the double peaks in the relevant 
thermodynamic response functions. 

A presence of the intermediate magnetization plateau 
has also been reported for the hybrid diamond-chain 
model, in which the spins of mobile electrons delocalized 
over decorating quantum dimers (interstitial sites) inter¬ 
act with each other as well as with the localized Ising 
spins situated at nodal lattice sites (Isl-IE^- The hybrid 
spin-electron diamond chain with localized Ising spins 
and delocalized electrons has been exactly solved using 
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an exact diagonalization procedure in combination with 
the decoration-iteration mapping technique. In this class 
of models, the competition emerges from the quantum- 
mechanical hopping of the interstitial electrons and is 
termed as the kinetically-driven spin frustration 

While the quantum entanglement has been rather ex¬ 
tensively studied in the Ising-Heisenberg diamond chains 
[HsMHBjl. it has been almost untouched in the analogous 
spin-electron diamond chains. It has been evidenced 
for the spin-1/2 Ising-Heisenberg diamond chain that 
the ground state is disentangled (classical) when the 
Ising coupling between the nodal and interstitial spins 
is the predominant one, but it becomes entangled for a 
strong enough Heisenberg coupling between the intersti¬ 
tial spins. The quantum concurrence generally decreases 
due to thermal fluctuations until it completely vanishes 
above a certain threshold temperature. The dependence 
of the quantum concurrence on the temperature and ex¬ 
ternal field was explored for example in the Ising-XXZ 
|5^ as well as Ising-XYZ [55| diamond chain. In the lat¬ 
ter model, the XY anisotropy may result in a re-entrant 
behavior of the quantum entanglement. 

Although thermodynamic properties of the hybrid 
spin-electron diamond chain have been already explored 
in some detail the inter-relation between the 

kinetically-driven spin frustration and the fermionic 
quantum entanglement is still missing [56|. In this work, 
we will address this question by considering an exactly 
solvable diamond chain with localized Ising spins and mo¬ 
bile electrons as a prototype model. The kinetic term 
associated with the electron mobility generates compe¬ 
tition and spin frustration. We will report the ground 
state phase diagram including ferrimagnetic (FRI), sat¬ 
urated paramagnetic (SPA) and unsaturated paramag¬ 
netic (UPA) phases. We will examine in detail the 
magnetization process, which exhibits an intermediate 
plateau, as well as the spin correlations around the pla- 
quette. A frustrated regime will be characterized by the 
absence of a local order that simultaneously satisfy all 
first-neighbors correlations. We will also compute the 
fermionic quantum concurrence between a pair of inter¬ 
stitial electrons in order to evaluate the influence of ther¬ 
mal fluctuations in the degree of quantum entanglement 
and its relation with spin frustration. We will unveil dis¬ 
tinct regimes of quantum entanglement as a function of 
the hopping amplitude and the external field, including 
the emergence of a re-entrant fermionic concurrence in 
the vicinity of the zero-temperature FRI-UPA and SPA- 
UPA ground-state phase transitions. 

This work is organized as follows: In section H, we will 
describe the model Hamiltonian and present the method¬ 
ology used to obtain the exact solution for the magneti¬ 
zation, correlation functions and quantum concurrence. 
In section H, we present the ground state phase diagram, 
and a detailed study of the magnetization processes, cor¬ 
relation functions, spin frustration and quantum entan¬ 
glement. In section IV, we summarize and draw our main 
conclusions. Some details of the analytical derivation are 


given in the appendices. 


II. MODEL AND ITS HAMILTONIAN 

Let us consider a coupled spin-electron model on a dia¬ 
mond chain, which contains localized Ising spins = ±1 
on its nodal sites and two mobile electrons on each couple 
of interstitial sites (see Fig. [1]). For further convenience, 
the total Hamiltonian can be defined as a sum over cell 
Hamiltonians, i.e. 1-L = whereas each cell Hamil¬ 

tonian Hi involves all the interaction terms belonging to 
Tth diamond unit: 

2 

'Hi = -t (4^ ,^Ci 2,7 + h.c.) 

7=t,i i=i 

2 

+ J (di -h CFi+l) ~ 

i=i 

- 2 (^) 

Here, ct^-^ and Cij^j {j = 1 , 2 ) are the usual fermionic 
creation and annihilation operators for the mobile elec¬ 
trons with the spin 7 or ^ and riij = is the 

respective number operator. The hopping term t takes 
into account the kinetic energy related to a quantum- 
mechanical hopping of the mobile electrons on the in¬ 
terstitial sites, the coupling constant J accounts for the 
nearest-neighbor Ising interaction between the localized 
Ising spins and the mobile electrons, and h is Zeeman’s 
energy of the localized Ising spins and mobile electrons 
in a presence of the external magnetic field. 


Itl> ltlp> ltt> 



FIG. 1: A diagrammatic representation of a correlated spin- 
electron model on a diamond chain. The nodal sites are occu¬ 
pied by the localized Ising spins ai = ±1, which are coupled 
by the Ising interaction with the spin of mobile electrons hop¬ 
ping between the nearest-neighbor interstitial sites. 


The matrix form of the Hamiltonian (HD in the local 
basis of two mobile electrons for the i-th diamond cell 
|4-4)i, Itr0)i, It4)i, |0,t4-)i and can be 
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represented as: 
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where we have introduced the notation for the total spin 
of two Ising spins from i-th diamond unit jii = +(7^+1. 

For simplicity, we have left out the constant term —h[iij2 
as it can be later simply added to the corresponding 
eigenvalues of the Hamiltonian matrix (j2j). A straight¬ 
forward diagonalization of the Hamiltonian matrix ([2]) 
leads (after taking into account also the constant term 
—hiii/2) to the following spectrum of eigenvalues: 


= 

2J III 2h 2 ^*’ 

( 3 ) 

^i2 = 

2JiJji -j- 2h 

( 4 ) 


h 

— “Mz + 2t, 

( 5 ) 
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( 6 ) 
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( 7 ) 
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( 8 ) 

— 



which correspond to the eigenvectors 

l¥’i2) 14-5 4')i ) 


( 9 ) 

( 10 ) 


Wi'i) — + lOitDi + ItiiO),], (11) 

Wii) = ^[ii,t)i + it4)i-|o,n)i-in,o)j, (12) 


1 ^, 5 ) = ^[|o,n),-in,o)j, 


(13) 

(14) 


All physical quantities of interest can be subsequently 
calculated from these exact analytical results by adapting 
the standard transfer-matrix method. 


A. Magnetization and correlation functions 


matrix technique. The partition function of the corre¬ 
lated spin-electron model on a diamond chain can be cal¬ 
culated following the procedure: 


N N 6 

{a}i=l {a}i=lj=l 

N 

{a} i=l 


where [3 = 1/(/cbT), ks is Boltzmann’s constant, T is 
the absolute temperature, the summation carried 

out over all possible spin configurations of the localized 
Ising spins, the symbol Tr^ refers to the trace over the 
degrees of freedom of two mobile electrons from the i-th 
diamond unit and W = The 

expression W can be regarded as the transfer matrix: 

^ (u:{2) u;(0) \ 

Vc^(O) uj{-2) )' 

whereas individual matrix elements depend just on the 
total spin of two localized Ising spins from i-th diamond 
unit according to: 



w(/ii) = w((Ti,(Ti+i) = ^exp(-/3eij) = 2exp /iij 

X [1 + cosh(2/3t) + cosh(2/3J/ii — 2ph)] . (17) 

The exact result for the partition function ([15]) then read¬ 
ily follows from the two eigenvalues of the transfer matrix 

(HD: 


A± — - [^(2) + ^(—2) ± Q] , (18) 

where Q = ([ci;(2) — ci;(—2)]^ -h 4[ci;(0)]^) ^. In the ther¬ 
modynamic limit of A ^ (X), only the largest eigenvalue 
effectively contributes to the partition function. 

The unitary transformation that diagonalizes the 
transfer matrix (HD, is determined by the matrices U 
and U~^\ 


/ A+-a;(-2) A_-u;(-2) \ 

U= ^ (19) 

\ A+ / 

where A± = y^[A± — uj{—2)Y + [^(0)]^ and U~^ = I/'*'. 

Now, one may take advantage of the unitary trans¬ 
formation in order to calculate the magnetization and 
correlation functions. For instance, the magnetization of 
localized Ising spins at the nodal lattice sites follows from 
the relation: 


Here, we will derive exact results for the magnetization 
and correlation functions, as provided by the transfer- 


Ma = {ct!) = ^Tr 




N 


( 20 ) 
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where = U~^a^U and W = U~^WU. is the stan¬ 
dard 2 x 2 Pauli matrix associated with the z-component 
of a nodal Ising spin. A similar formula can also be de¬ 
rived for the magnetization at the interstitial site j = 1,2 
of the i-th cell: 


Ms = {S!j) = —Tr 




( 21 ) 


where = U with the elements (T,+i) = 

Tr and the trace performed over the states of 

the interstitial sites belonging to the i-th cell. The matrix 
S^j corresponds to the z-component of the spin operator 
at the interstitial site j = 1 , 2 of the i-th cell, written in 
the corresponding 6 x 6 sub-space. The total magnetiza¬ 
tion can be then obtained either by adding both individ¬ 
ual contributions of the localized nodal Ising spins (|2Q|) 
and the mobile electrons (ED, or respectively, from the 
usual formula: 


M = 


(3 dh 


InZ 


AT- 


( 22 ) 


The same procedure can be used to calculate the correla¬ 
tion functions. The correlation functions between various 
spatial components of spin operators S'g- at the intersti¬ 
tial sites of an elementary plaquette can be calculated as 


= ^Tr , (23) 

where with the elements 

^i+i) = Tr Sf-^Sf 2 e~^^K The correlation 

function between the localized Ising spins and the 
interstitial spins is given by: 






(24) 


B. Fermionic entanglement 


Quantum entanglement is closely related to non-local 
correlations present in a quantum system. For two 
qubits, it can be quantified by the entanglement of for¬ 
mation 1 ^: 


where 


Ef=H 




(25) 


H{x) = -xXog^ix) - (1 - a:;) log 2 (l - x), (26) 

and C is called the quantum concurrence given by: 

C = maxjo, \/a 7-vT-- V^l • (27) 

Here, are the eigenvalues of the matrix R = 
p(cr^ (g) (jy) p* (cr^ (g) (7^) sorted in descending order, 


is the usual Pauli matrix and p represents the reduced 
density matrix for a pair of qubits. Because the entan¬ 
glement of formation (j25j) is a monotonous function of the 
concurrence C, one may directly use the concurrence as 
a measure of quantum entanglement ranging from 0 (no 
entanglement) up to I (maximum entanglement) 

In the present model, we will quantify the quantum 
entanglement in the sub-space of up spins in a given pla¬ 
quette [m. Therefore, in the local basis |0, 0), |0, t), It: 0) 
and |t:t): fhe reduced density matrix can be written in 
the general form 


/ a 0 0 0 \ 

Ox z 0 I 
0 z* p 0 ’ 

Vo 0 06/ 


(28) 


whose elements can be directly associated with the fol¬ 
lowing thermal averages 

b = (nii,t,«i 2 ,t) (29) 

X = (na,t) - 

y = (na,t) - (^ii,t,^i 2 ,t) 

a = ((1 - nii.t)(l - ni 2 ,t)) 

= V = 

Following the procedure elaborated previously in Ref. 
[m, one arrives to the following explicit form of the 
eigenvalues: 

Al = (|(nil,t) - (nil.tni2,t)l + K4l,Vi2,t)l) , 

A 2 = (|(nii,t) - (nii,tni 2 ,t)l - l(4i,Vi2,t)l) , 

A 3 ,4 = (na,t«i 2 ,t) (1 - 2(na,t) + {nii,i-ni2,t)) ■ (30) 

According to Eqs. pip . the fermionic concurrence can 
be calculated from the relation: 

C = 2 maxjo, 1(4 ,,,Cj 2 ,t)| 

- /(«ii,t^i2,t) (1 - 2(nii,t) + (nii.tni2,t))| (31) 

which depends on the statistical mean values 

and explicitly given in Appendix B. It 

is noteworthy that the more general formula (|3T|) for the 
fermionic concurrence reduces to the simpler formula re¬ 
ported previously by Deng and Gu [58| when considering 
the particular case {nn^^) = (^^ 24 ) = 1 / 2 . 


III. RESULTS AND DISCUSSIONS 

Now, let us proceed to a discussion of the most in¬ 
teresting results for the coupled spin-electron diamond 
chain. For simplicity, our further discussion will be re¬ 
stricted to a particular case with the antiferromagnetic 
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Ising interaction J > 0, which should exhibit the most 
diverse magnetic behavior on behalf of a mutual inter¬ 
play between the kinetically-driven spin frustration, local 
quantum fluctuations and the Zeeman’s splitting. 

A. Ground state 

First, we will turn our attention to the ground-state 
phase diagram shown in Fig. [21 which involves the ferri- 
magnetic (FRI), the saturated paramagnetic (SPA) and 
the unsaturated paramagnetic (UPA) ground states un¬ 
ambiguously determined by the eigenvectors: 

N 

iFRi) = ni^ii)®iai=-i), 

i=l 

N 

|SPA) = ® |(Ti = 1), 

i=l 

N 

|UPA) = n 1^*4) ® ki = 1), (32) 

i=l 

with the corresponding eigenenergies per plaquette: 

^FRi = -4 J - h, 

^SPA = 4 J - 3h, 

F’upa = —h — 2t. (33) 

It can be understood from the eigenvectors (|3^ that the 
interstitial spins are aligned into the magnetic field and 
the nodal spins in opposite to the magnetic field within 
the FRI ground state, while all nodal as well as inter¬ 
stitial spins are equally aligned into the magnetic field 
within the SPA ground state. However, the most peculiar 
spin arrangement can be found within the UPA ground 
state, where the nodal spins are aligned into the mag¬ 
netic field but the interstitial spins are subject to a quan¬ 
tum entanglement of two antiferromagnetic and two ionic 
states as described by the eigenvector m- The phase 
UPA originates from a kinetically-driven spin frustration, 
which is closely connected to the quantum-mechanical 
hopping of two mobile electrons with opposite spins on 
the interstitial sites. In agreement with this statement, 
the nodal spins are completely free to flip in an absence 
of the magnetic field and the UPA ground state becomes 
highly degenerate. High macroscopic degeneracy can be 
also found at a triple point given by the coordinates 
t/J = 2.0 and h/J = 4.0, at which all available ground 
states FRI, SPA and UPA coexist together with another 
special nodal antiferromagnetic (NAF) ground state (see 
Appendix A). 

B. Correlation functions and magnetization 

The local spin arrangements inherent to each individ¬ 
ual ground state can be witnessed through the pair cor¬ 
relation functions, which can also give insights into how 



FIG. 2: The ground-state phase diagram in the t/J vs. h/J 
plane. In the UPA phase, the interstitial electrons are quan¬ 
tum entangled. It is highly degenerated at zero field because 
the nodal spins become effectively decoupled. A high degen¬ 
eracy is also present at the triple point with the coordinates 
t/J — 2.0 and h/J — 4.0, where all three phases co-exist with 
a nodal antiferromagnetic phase (see Appendix A). 

the spin-spin correlations are affected by thermal fluctu¬ 
ations. 



kgT/J kgT/J 


FIG. 3: Temperature variations of the transverse correla¬ 
tion between the interstitial spins {S /{^the longitudinal 
correlation between the interstitial spins {S/iS/ 2 )^ the longi¬ 
tudinal correlation between the interstitial and nodal spins 
(crhS'fi) for the set of parameters: (a) t/J — 1.5, h/J — 3.0 
(FRI phase); (b) t/J = 2.5, h/J = 3.0 (UPA phase); (c) 
t/J = 1.5, h/J = 5.0 (SPA phase); (d) t/J = 2.0, h/J = 4.0 
(triple point). Notice that spin frustration is present in the 
low temperature regime above the UPA ground state where 
only out of three longitudinal spin-spin correlations within 
the elementary plaquette is negative. 

In Fig. [3] we have plotted, for all available ground 
states, typical thermal dependences of the transverse cor¬ 
relation function between the interstitial spins {SfiSf 2 ) = 
(S'fiS'f 2 ) (dash-dot lines), the longitudinal correlation 
function between the interstitial spins {S/ 1 S/ 2 ) (dash 
lines) and the longitudinal correlation function between 
the nodal and interstitial spins (dash-dot-dot 

lines). The zero-temperature asymptotic limits of the 
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correlation functions shown in Fig. [Sfa) clearly con¬ 
firm anti-parallel spin alignment between the interstitial 
and nodal spins and the parallel alignment of intersti¬ 
tial spins within the classical FRI phase, whereas an in¬ 
crease in temperature may be responsible for an up-rise 
of small transverse (quantum) correlations. It is clear 
from Fig. [Sj^b) that the longitudinal correlation between 
the interstitial spins {SfiSf 2 ) is antiferromagnetic at low 
but ferromagnetic at high temperatures, which indicates 
a spin frustration above the UFA ground state accord¬ 
ing to the concept of temperature-dependent spin frus¬ 
tration [59|. The spin correlations within the interstitial 
spins do not reach unity due to the presence of strong 
quantum fiuctuations associated with the hopping pro¬ 
cess. The complete spin alignment of the interstitial 
and nodal spins relevant to the SPA ground state is evi¬ 
dent from the correlation functions depicted in Fig.ISfc), 
where a small increase in transverse correlation function 
induced by thermal fiuctuations can be also observed. 
Finally, Fig. ^d) illustrates the correlation functions 
exactly at the highly degenerate triple point, at which 
zero temperature asymptotic limits reach nontrivial val- 
ues = (5 - V5)/20, = (5 + 3^5)/20 

and (crfiSfi) = (3\/5 — 5)/10 (see Appendix A). 



k„T/J k„T/J 

a D 


FIG. 4: Temperature dependence of the single-site magneti¬ 
zation Ms of the interstitial spins, the single-site magnetiza¬ 
tion Ma of the nodal spins and the total magnetization (M) 
per diamond unit for the set of parameters: (a) tjJ— 1.5, 
h/J = 3.0 (FRI phase); (b) t/J = 2.5, h/J = 3.0 (UPA 
phase); (c) t/J = 1.5, h/J — 5.0 (SPA phase); (d) t/J — 2.0, 
h/J = 4.0 (triple point). 


Thermal behavior of the total magnetization M per di¬ 
amond unit (dash-dot-dot lines) is plotted in Fig.|4]along 
with the single-site magnetization Ms of the interstitial 
spins (dash-dot lines) and the single-site magnetization 
Mfj of the nodal spins (dash lines). Fig. SJ^a) serves in 
evidence that the rising temperature may facilitate a spin 
reversal of the nodal spins within the FRI phase, which 
in turn causes a unusual thermally-induced increase in 


the total magnetization. Thermal excitations from the 
UPA ground state may also give rise to a striking in¬ 
crease of the total magnetization, which occurs due to 
a spin reorientation of the interstitial rather than the 
nodal spins (see Fig. [4|(b)). Contrary to this, the mag¬ 
netization of the nodal as well as interstitial spins show 
only a smooth monotonous decrease with the increas¬ 
ing temperature when starting from the SPA ground 
state (see Fig. S^c)). Last but not least, the magneti¬ 
zation curves for the highly degenerate triple point is 
illustrated in Fig. SJ^d). These start from the nontrivial 
zero-temperature asymptotic values Ms = (5 + \/5)/10, 
Mcj = v^/5 and M = (5 + 2^/b)/b (see Appendix A). 



FIG. 5: The field dependence of the magnetization for a few 
different values of the temperature and two selected values of 
the hopping term t/| J| = 1.5 and 2.5. The upper panel shows 
the single-site magnetization of the nodal spins, the central 
panel the single-site magnetization of the interstitial spins and 
the lower panel the total magnetization per diamond unit. 

Both available magnetization scenarios are displayed in 
Fig. [5l The magnetization curves shown in Figs. [5fa)- 
(c) shed light on a field-induced transition from the FRI 
phase to the SPA state, which occurs just for a sufficiently 
small values of the hopping term t/J < 2. Strictly speak¬ 
ing, the true discontinuous phase transition associated 
with the magnetization jump of the nodal spins occurs 
at zero temperature only (see Fig. [Sfa)). However, it 
is worthwhile to remark that the tendency of the nodal 
spins to align in opposite to the magnetic field is clearly 
visible at non-zero temperatures as well. A sudden reori¬ 
entation of the nodal spins in a vicinity of the saturation 
field causes a subtle decline of the magnetization of the 
interstitial spins due to the antiferromagnetic coupling 
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between the nodal and interstitial spins (see the inset 
in Fig. Efb)). Another magnetization scenario shown 
in Figs. [5fd)-(f) relates to a discontinuous field-induced 
transition from the UFA state to the SPA state. It can 
be clearly seen from Fig. [5fe) that the observed magne¬ 
tization jump near the saturation field occurs on account 
of a spin reversal of the interstitial spins. The sudden 
change in the magnetization of the interstitial spins now 
evokes a subtle decline of the magnetization of the nodal 
spins. Altogether, it could be concluded that the mag¬ 
netization curve always displays a magnetization plateau 
at one-third of the saturation magnetization even though 
the microscopic mechanism for a plateau formation may 
be different. 


C. Spin frustration 


It is quite apparent from the eigenvector (|3^ of the 
UPA ground state that the quantum-mechanical motion 
of two mobile electrons on the interstitial sites gives rise 
to a kinetically-induced spin frustration of the nodal Ising 
spins at zero temperature. To verify a frustrated charac¬ 
ter of the nodal spins at finite temperatures one may take 
advantage of the concept of temperature-dependent frus¬ 
tration, which requires a negative sign for the product of 
correlation functions along an elementary plaquette im¬ 
plying incapability of the spins to satisfy all underlying 
spin-spin interactions. The negative sign of the prod¬ 
uct {S^iS^ 2 ) t)e achieved just if the longitu¬ 

dinal correlation function between the interstitial spins 
{SkSh) is predominantly antiferromagnetic (negative). 
It has been demonstrated in the previous section that 
the correlation function indeed changes its sign 

at a certain temperature when starting from the UPA 
ground state, whereas the relevant frustration temper¬ 
ature can be regarded as an indicator of the transition 
from the frustrated regime to the non-frustrated one. 

The frustration temperature, as calculated from cross¬ 
ing points of the longitudinal correlation function be¬ 
tween the interstitial spins is plotted in Fig. [6] 

against the hopping term t/J for a few selected values 
of the magnetic field. It turns out that the antiferro¬ 
magnetic correlation { 8 ^ 18 ^ 2 ) ^ ^ develops just for a 
sufficiently strong hopping term t/J > 2, because the 
coupled spin-electron diamond chain may become frus¬ 
trated only within the UPA phase. The frustration tem¬ 
perature accordingly starts at low enough magnetic fields 
h/J<A from the ground-state boundary between the 
FRI and UPA phases at t/J = 2, while the initial value 
of the frustration temperature is progressively shifted 
along the ground-state boundary between the SPA and 
UPA phases t/J = hjJ — 2 at higher magnetic fields 
h/J > 4. Altogether, it can be concluded that the frus¬ 
trated regime of the correlated spin-electron diamond 
chain is delimited from above by the displayed curves 
of frustration temperature. 



FIG. 6: The frustration temperature as a function of the hop¬ 
ping term for a few selected values of the magnetic field. The 
investigated system is frustrated below the displayed lines 
where the product of the longitudinal spin correlations around 
a plaquette is negative. Frustration is only present above the 
UPA ground state. 


D. Fermionic concurrence 


Last but not least, we will turn our attention to a dis¬ 
cussion of the bipartite quantum entanglement between 
two mobile electrons located at the interstitial sites of a 
given plaquette, which will be quantified by the fermionic 
concurrence calculated according to Eq. dSID. The tem¬ 
perature dependence of the fermionic concurrence is de¬ 
picted in Fig.[7]for several values of the hopping term and 
two different values of the magnetic field. As it could be 
expected, the concurrence is zero well within the classical 
FRI phase as illustrated by the solid line corresponding to 
t/J = 1.0 in Fig. [71^a). On the other hand, the fermionic 
concurrence within the fully entangled UPA ground state 
for t/J = 3.0 (dash-dot-dot line) starts from its maxi¬ 
mum value and then gradually decreases with increasing 
temperatures. The special case with the fixed value of 
the hopping term t/ J = 2.0 corresponds to the phase co¬ 
existence of the FRI and UPA phases at zero temperature 
(the dashed line) and hence, the zero- temperature limit 
of the fermionic concurrence is given by the mean value of 
the disentangled FRI phase and the fully entangled UPA 
phase. However, the most striking finding concerns with 
a re-entrant behavior of the fermionic concurrence when 
the hopping parameter is sufficiently close but slightly 
below the ground-state boundary between the FRI and 
UPA phases. Under this condition, the fermionic con¬ 
currence starts from zero in agreement with the classical 
character of the FRI ground state, then it emerges at 
some lower threshold temperature due to thermal excita¬ 
tions to the UPA phase and finally, it completely disap¬ 
pears at a certain higher threshold temperature (see the 
dash-dot line for the fixed value tjJ = 1.5) . 

It can be understood from Fig.[7|(b) that the fermionic 
concurrence displays qualitatively similar thermal varia- 
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FIG. 7: The fermionic concurrence versus temperature for 
several values of the hopping term and two different values of 
the magnetic field: (a) hjJ — 3.0 (transition from the FRI 
to the UFA ground state); (b) h/J — 5.0 (transition from 
the SPA to the UPA ground state). The re-entrant behavior 
emerges when the hopping amplitude drives the investigated 
system to the classical FRI or SPA ground state, but still 
keeps it in a close vicinity of the ground-state boundary with 
the quantum UPA ground state. 

tions also at higher magnetic fields when the variation of 
the hopping term causes a phase transition between the 
UPA and SPA ground states. As one can see, the concur¬ 
rence always starts from zero in accordance with the clas¬ 
sical character of the SPA ground state for small enough 
hopping terms, but it may eventually show the interest¬ 
ing re-entrant behavior when the hopping term is selected 
sufficiently close to the ground-state boundary between 
the SPA and UPA phases (see dash-dot line for the hop¬ 
ping term tjJ = 2.0 in Fig.[71^b)). The dashed line shows 
the fermionic concurrence exactly at the ground-state 
boundary between the SPA and UPA phases (t/J = 3.0 
ioT hjJ = 5.0) when the concurrence starts from its mean 
value C = 1/2. If the hopping term is strong enough to 
enforce the UPA ground state, then, the fermionic con¬ 
currence starts from its maximum value and afterward 
gradually diminishes with the rising of temperature until 
it completely vanishes above a certain threshold temper¬ 
ature. 

The threshold temperature for the disappearance of 
the fermionic concurrence is plotted against the hopping 
term in Fig. [8] for three different values of the magnetic 
field. The mobile electrons on the interstitial sites are 
fully entangled just within the UPA phase and hence, it 
follows that all depicted lines of threshold temperature 
should start from the ground state boundary between the 
UPA-FRI and UPA-SPA phases, respectively. The de¬ 
pendences presented in Fig. [8] are indeed consistent with 
this statement and they also provide an independent con¬ 
firmation of the striking re-entrance when the quantum 
entanglement is induced above the classical FRI or SPA 
phase on account of vigorous thermal excitations to the 
UPA state. 

At this stage, it might be quite interesting to study 
a mutual correlation between the spin frustration and 
quantum entanglement, because the kinetically-driven 
spin frustration of the nodal spins appears due to the 



FIG. 8: The threshold temperature as a function of the hop¬ 
ping amplitude for three different values of the magnetic field: 
h/J = 0.0, 3.0 and 5.0. The system is entangled (disen¬ 
tangled) inside (outside) the area delimited by the displayed 
lines. The re-entrant behavior indicates that thermal fluctu¬ 
ations favors quantum entanglement in the close vicinity of 
the FRI-UPA and SPA-UPA ground-state phase boundaries. 




FIG. 9: A comparison between the threshold and frustra¬ 
tion temperatures plotted as a function of the hopping am¬ 
plitude for two different magnetic fields: (a) h/J = 0.0; 
(b) h/J — 5.0. The notation for parameter regions are: 
Entangled (E), Erustrated (E), Disentangled (D), and Non- 
frustrated (NE). The re-entrance in the threshold tempera¬ 
ture depicted in Eig. [9ja) implies the emergence of a low- 
temperature entangled-non-frustrated regime, as well as, a 
high-temperature disentangled-frustrated regime. 


hopping process of the mobile electrons with opposite 
spins inevitably underlying the quantum entanglement. 
Bearing this in mind, the fermionic quantum entangle¬ 
ment could be thus regarded as an indispensable ground 
for the kinetically-driven frustration of the nodal spins. 
It actually turns out that the zero-temperature asymp¬ 
totic value of the threshold and frustration temperature 
always start from the same value of the hopping term, 
whereas the threshold temperature for the quantum en¬ 
tanglement generally lies slightly above the frustration 
temperature at low and moderate temperatures. How¬ 
ever, the threshold and frustration temperature might 
also cross each other at a certain higher temperature as il¬ 
lustrated in Fig. [9ja), which gives evidence that the spin 
frustration of nodal spins may emerge at higher temper- 
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ature despite of disentangled character of the interstitial 
spins. However, we must have in mind that the crossing 
of the frustration and threshold temperature takes place 
at a certain temperature, which is however high enough 
to produce only very weak spin-spin correlations around 
the elementary plaquette. 


IV. SUMMARY AND CONCLUSIONS 

In the present work, we have provided a detailed analy¬ 
sis of the relationship between spin frustration and quan¬ 
tum entanglement in a hybrid spin-electron system con¬ 
taining the localized Ising spins and mobile electrons. In 
particular, we have considered the diamond chain proto¬ 
type model with the localized Ising spins situated at the 
nodal lattice sites that interact with two delocalized elec¬ 
trons quantum-mechanically hopping between a pair of 
the interstitial sites. The model was exactly solved using 
the transfer-matrix technique. The ground-state phase 
diagram was shown to be composed of two classical SPA 
and FRI phases in addition to the one quantum UPA 
phase, in which the interstitial electrons are quantum- 
mechanically entangled. The relevant spin-spin correla¬ 
tions around an elementary plaquette were obtained and 
unveiled a low-temperature frustrated regime in the pa¬ 
rameter region corresponding to the UPA ground state. 
In this regime there is no local spin ordering that satisfies 
simultaneously all signs of the spin correlations. In the 
present model frustration results from the quantum hop¬ 
ping process that favors the anti-parallel spin alignment 
of the interstitial electrons. We have also analyzed the 
distinct magnetization processes and showed the occur¬ 
rence of intermediate plateau at one-third of the satura¬ 
tion magnetization. Also, we have found that the spin re¬ 
orientation of one sub-lattice (nodal or interstitial spins) 
observable in a vicinity of the field-induced transition is 
responsible for a small drop in the magnetization of the 
complimentary sub-lattice. Finally, we have quantified 
the degree of fermionic entanglement between two inter¬ 
stitial electrons by computing the quantum concurrence 
in the spin up sub-space. 

The low-temperature asymptotic value of the fermionic 
concurrence has unambiguously confirmed an existence of 
the quantum entanglement within the UPA ground state, 
while it has also demonstrated rather unexpected re¬ 
entrant behavior due to thermally activated excitations 
when the hopping term drives the investigated model 
towards the disentangled FRI and SPA ground states 
but still preserves it in a vicinity of the ground-state 
boundary with the quantum UPA ground state. We have 
also compared the threshold and frustration tempera¬ 
tures with the goal to bring insight into a mutual rela¬ 
tionship between the quantum entanglement and geomet¬ 
ric spin frustration. It has been convincingly evidenced 
that these two physical properties are strongly related 
yet independent. Besides the conventional entangled- 
frustrated and disentangled-non-frustrated regimes, we 


have additionally found entangled-non-frustrated and 
disentangled-frustrated regimes as well. Although geo¬ 
metric spin frustration and quantum entanglement are 
closely related features within the present model, this 
result indicates that thermal fluctuations have distinct 
influence on these two physical aspects. It would be valu¬ 
able to have further studies of quantum spin chains with 
competing interactions to shed additional light of the rel¬ 
evant thermal processes influencing the interplay between 
spin frustration and quantum entanglement. 
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VI. APPENDIX A 



FIG. 10: Degenerated plaquette configurations at the co¬ 
existence point with the coordinates /i/J = 4 and tjJ — 2: 
(a) daturated paramagnetic state (SPA) with the interstitial 
and nodal spins parallel to the external field; (b) unsaturated 
paramagnetic state (UPA) with nodal spins parallel to the ex¬ 
ternal field and interstitial electrons antiparallel to each other 
in an entangled state; (c) ferrimagnetic state (FRI) with the 
nodal spins anti-parallel to the field and interstitial spins par¬ 
allel to the external field; (d) nodal antiferromagnetic state 
(NAF) with the nodal spins anti-parallel with respect to each 
other and interstitial spins parallel to the external field. 


At the point of co-existence between the SPA, FRI and 
UPA phases, their corresponding configurations in a pla¬ 
quette are also degenerated with a nodal antiferromag¬ 
netic phase (NAF). Therefore, in the ground state, cells 
with either one of these configurations are allowed (see 
Fig. 10). However, NAF cells can only appear in pairs 
and only FRI cells can be placed between them. Defining 
Nspa^ NupAi and N^af as the number of cells in 

each configuration {N = Nspa-\-Nupa-\-Nfri^Nnaf) 
the number of possible ways to distribute then along the 
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chain can be written as 

^ _ {NsPA + Nupa)\ ^ 

Nspa^‘Nupa\ 

{Nspa + Nupa + Nnaf/^NjI ^ 

{Nspa + Nupa)1 {Nnaf/^NjI 

{Nnaf/‘^ + Nfri)\ 

{NNAF/‘Nj\NFRi\ 

where the first fraction accounts for the distinct permu¬ 
tations between SPA and UP A cells, the second fraction 
accounts for the number of ways NAF dimers and para¬ 
magnetic cells (either SPA or UPA) can be arranged, and 
the last one corresponds to the number of possible ways 
FRI cells can be distributed between NAF dimers. The 
actual fraction of cells in each configuration is found by 
maximizing fl, resulting in 


Nspa/N = 

5 - x/5 

10 ’ 

(35) 

Nupa/N = 

5- \/5 

10 ’ 

(36) 

Nfri/N = 

10-4\/5 

10 

(37) 

Nnaf/N = 

-10 + 6\/5 

10 

(38) 


The values of magnetization of the nodal and interstitial 
sites, as well as the relevant correlation functions in each 
configuration are summarized in table L The proper aver¬ 
age of these quantities, weighted considering the relative 
fraction of cells in each configuration, provides the values 
reported in the main text at the triple co-existence point 
h/J = 4 and t/J = 2. 


a function of the proper reduced matrix. One starts by 
defining the auxiliar operator associated with the eigen¬ 
states of a pair of interstitial electrons : 


6 

i=i 


which depends on the spin orientation of the neighboring 
Ising spins. The reduced density matrix elements can be 
obtained by tracing out all Ising spins and all interstitial 
electrons except those at a given plaquette. It can be 
directly shown that the reduced density matrix elements 
can be put into the form: 



y](w ((71,(72).. 


{o'r ~\~ (7r+l) 


X w((7r+l,(7r+2) • • . w((7iv, (7i)) 


(40) 


where we introduced the matrix 


P 


k,l = 


^/c,z(2) gk,i{0) A 

Qk,l{^) Qk,l{-‘^) J ' 


(41) 


and the indices k and I run from 1 to 6 according to 
the possible states of the interstitial sites 
Itr0)i, It4)i, |0,t4.)i and |4,,t)i. 

By exploring the unitary transformation that diago¬ 
nalizes the transfer matrix uj (see main text), the density 
matrix elements result in 


TABLE I: Fraction of cells Nx/N (X=SPA, UPA, FRI, NAF) 
in each cell configuration at the co-existence point {hjJ — 4, 
t/J = 2). The corresponding values of the nodal magnetiza¬ 
tion M(j^ interstitial magnetization Ms, correlation between 
nodal and interstitial spins {(JiSfi), correlation between the 
spin ^-components of the interstitial sites {Sfi sh), and be¬ 
tween their spin x-components (AS'fiAS'f 2 ) are also shown. Their 
weighted average provides the values reported in the main 
text. 


Plaquette 

Nx/N 


Ms 


^ {ShS!^) 

(SfiSf,) 

SPA 

5-V5 

10 

1 

1 

1 

1 

0 

UPA 

5-\/5 

10 

1 

0 

0 

- 1/2 

1/2 

FRI 

IO- 4 U 5 

10 

-1 

1 

-1 

1 

0 

NAF 

-IO+ 6 U 5 

10 

0 

1 

0 

1 

0 


Pk,l = 


Qk,i{‘^) + Qk,i{—‘^) 2pfe^/(Q)u;(Q) 

2A+ QX-\- 

, [eM(2)-eM(-2)][^(2)-c^(-2)] 


2QA+ 


(42) 


As a function of the reduced density matrix elements, 
one can readly obtain the thermodynamic averages 

= 2/)5,4 (43) 

= Pl,l 

= Pl,l + 2/93^3 

where the identities ^ 5^4 = p 6,3 and = Pa,a were used. 
From the above expressions, the quantum concurrence 
can be explicitly written as 


VII. APPENDIX B , _. 

C = max I 0,4/95,4 - 2^/0i,i(l - /9i,i - 4 / 93 , 3 ) t (44) 

The relevant thermal averages of quantities associated ^ ^ 

with a pair of interstitial electrons can be obtained as 
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